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ABSTRACT 

An  equation  of  the  K-dV  type  is  derived  for  the  study  of  nonlinear 
surface  waves  under  gravity  with  surface  tension  and  edge  constraints  in  a 
channel  of  variable  cross  section.  The  approach  used  here  is  based  upon  a 
nonlinear  ray  method.  An  example  of  a  symmetric  rectangular  channel  with 
variable  width  is  given. 
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SIGNIFICANCE  AND  EXPLANATION 

An  equation  of  the  K-dV  type  is  derived  for  the  study  of  nonlinear 
surface  waves  under  gravity  with  surface  tension  and  edge  constraints  in  a 
channel  of  variable  cross  section.  The  approach  used  here  is  based  upon  a 
nonlinear  ray  method.  An  example  of  a  symmetric  rectangular  channel  with 
variable  width  is  given. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


NONLINEAR  CAPILLARY  WAVES  UNDER  GRAVITY 
WITH  EDGE  CONSTRAINTS  IN  A  CHANNEL 


N.  C.  Shan* 

{ 1 .  INTRODUCTION 

The  dynamical  problem  of  capillary  waves  on  a  liquid  in  a  channel  has  not  been  much 
studied  in  the  past.  The  difficulty  may  lie  in  the  fact  that  the  edge  condition  at  the 
line  of  contact  of  the  liquid  free  surface  with  the  channel  wall  so  far  remains  an  open 
problem.  It  is  well  known  that  for  water  at  rest  the  contact  angle  at  the  edge  of  the  free 
surface  is  constant.  However/  when  the  surface  is  in  motion,  the  edge  condition  depends 
upon  whether  the  channel  wall  is  wetted  or  not  and  has  to  be  deterained  as  a  part  of  the 
solution  to  the  free  surface  problem.  In  our  study  of  capillary  waves  under  gravity  on  a 
viscous  fluid  in  an  inclined  channel  (Shih,  1973;  Shen  and  Shih,  1974),  we  made  the 
assumption  that  the  fluid  velocity  should  be  continuous  at  the  line  of  contact.  It  follows 
that  the  free  surface  elevation  will  remain  aero  if  initially  it  is  sero  there.  This 
argusMnt  would  seam  to  be  plausible  only  for  small  amplitude  waves.  By  a  different 
consideration  the  same  edge  condition  was  also  used  by  Kopachevskii  (1967). 

Zn  two  recent  papers  by  Benjamin  and  Scott  (1979)  and  Benjamin  (1981),  the  problem  of 
gravity-capillary  waves  with  edge  constraints  on  an  inviscid  fluid  has  been  investigated 
both  experimentally  and  theoretically.  They  considered  a  straight  horizontal  channel  of 
uniform  cross  section,  and  filled  the  channel  with  water  up  to  the  brim  so  that  the  edge 
condition  of  sero  surface  elevation  can  be  experimentally  realized.  In  the  former  paper, 
the  linear  problem  was  studied.  In  the  latter,  variational  principles  were  formulated  for 
the  full  nonlinear  equations  as  a  Hamiltonian  system  and  for  waves  of  permanent  type. 
However,  as  pointed  out  by  Benjamin  (1981),  a  rigorous  treatment  of  the  nonlinear  problem 
is  beyond  reach  at  present.  In  this  paper,  we  shall  use  a  different  approach  to  study  a 
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more  general  problem,  that  la,  the  channel  is  allowed  to  have  a  slowly  varying  cross 
section.  The  main  purpose  of  our  work  is  to  derive  an  equation  of  the  K-dV  type,  which  may 
be  used  as  an  approximate  equation  for  the  investigation  of  long  capillary  waves  with  edge 
constraints  in  a  channel. 

Our  approach  is  based  upon  the  nonlinear  ray  method  developed  in  Shen  and  Keller 
(1973).  The  basic  ideas  may  be  explained  as  follows.  The  motion  of  a  wave  front  is 
assumed  to  be  determined  by  a  phase  function.  It  is  found  that  the  phase  function 
satisfies  the  Hamilton-Jacobi  equation  in  geometrical  optics,  the  solutions  of  which 
determine  a  family  of  bicharacteristics  called  rays.  An  equation  of  the  K-dV  type  is  then 

derived  along  the  rays  for  the  wave  amplitude.  In  a  report  by  Zhong  and  Shen  (1982),  we 

specialized  the  procedure  in  Shen  and  Keller  (1973)  to  the  case  of  an  incompressible 
inviscid  fluid  without  surface  tension,  and  an  aquation  of  the  K-d V  equation  for  a  channel 
of  variable  cross  section  was  derived.  The  present  work  is  an  extension  of  the  previous 
result  and  a  host  of  problems,  such  as  fission  of  solitons  (Johnson,  1973»  Zhong  and  shen, 
19S2 ) ,  shelf  generation  behind  a  solitary  wave  (Knickbocker  and  Newell,  1980)  and  others, 
may  be  studied  similarly  by  means  of  the  equations  derived  here.  However,  this  report  is 
self-contained,  there  is  minimum  reference  to  the  previous  report.  He  also  note  that  our 
method  is  related  to  those  developed  in  Kuzmak  (1959)  and  choquet-Bruhat  (1969). 

He  formulate  the  problem  in  $2.  In  $3  we  derive  an  aquation  of  the  K-dV  type  for  a 

general  channel.  In  {4  the  example  of  a  symmetric  rectangular  channel  with  variable  width 
is  considered  to  illustrate  the  method.  Finally  some  discussions  are  given  in  |5. 


J2.  FORMULATION  OF  THE  PROBLEM 

H«  consider  the  notion  of  an  inviscid,  incompressible  fluid  of  constant  density  with 
surface  tension  in  a  channel  with  a  saooth  convex  boundary  defined  by  h*(x*,y*,r*)  -  0, 
where  t*  is  positive  upward  (Figure  1 ) •  The  lines  of  contact  are  two  curves  of 
intersection  of  h*  ■  0  with  the  plane  s*  •  d* .  The  governing  equations  are 

tti  *  VJ«  +  “  0  *  M) 

“t*  +  u*u5*  +  v*uy*  +  W*UJ*  "  -PJ./0*  *  (2) 

v£.  ♦  UJ.  ♦  v*vj,  ♦  WvJ,  -  -pJ./P*  .  (3) 

+  u<wj«  +  ♦  ***;•  •  -p;./p*  -  9*  ,  (4) 

subject  to  the  boundary  conditional  At 


s*  *  n*(x#,y*,t#>  , 

"t*  +  u*nJ*  +  n*.  -  w*  -0  , 


(5) 


u*hj#  ♦  v»h*.  ♦  w*h»#  -  0  .  (8) 

Here  (u*,v*,w*)  is  the  velocity,  t»  is  the  time,  g*  is  the  constant  gravitational 
acceleration,  P*  is  the  constant  density,  p*  is  the  pressure,  t*  m  n*  is  the  equation 


of  the  free  surface  and  T*  is  the  constant  surfaet  tension  coefficient.  Within  the 


framework  of  long  wave  approximation,  we  assume  that  the  channel  bottom  varies  slowly  in 
the  longitudinal  direction  and  the  magnitude  of  the  transverse  velocities  is  much  smaller 
than  that  of  the  longitudinal  velocity.  Under  these  assumptions,  we  introduce  the 
nondimens ional  variables 


-J- 


t  -  8-3/2tV(H/g)1/2,  (x,y,z)  -  (fl_3/2x*/H,g*/H,*VH)  . 

n  -  H*/H,  h  -  h*/H,  P  -  p*/(PgH)  . 

<u,v,e)  -  <u*/(gH)1/2,  B1/2v*/(gH>1/2,  B1/2w*/(gH) 1/2)  , 

T  -  T*/(PgH2),  B3/2  -  Vh  »  1,  d  -  d*/H  , 

where  L(H  are  raapactivaly  the  longitudinal  and  transverse  length  scales. 

the  unstarred  variables,  (1)  to  (8)  become 

u  +  6{v  +  w  )  -  0  , 
x  y  * 

ufc  ♦  uux  +  8(vuy  +  wu^)  -  -px  , 

v  +  uV  +  B(vu  +  vu  )  *  -B2p  , 
t  x  y  *  y 

w.  ♦  uw  +  8(vw  +  Ww  )  »  -B2(p  +  1)  , 
t  x  y  s  s 

subject  to  the  boundary  conditions: 

At  *  -  H(x,y,t)  , 

h.  +  un  +  B{vn  -  w)  ■  o  , 
w  x  y 

p  -  -TlfT^d  +  „*>  Mn(1  +  B“V) 
-28'3nxynxny] ( i  +  B"V  +  n2)-3/2  » 

at  h(x,y,c,s  -  d)  -  0  , 

h  -  d  ; 


In  teras  of 


(9) 

(10) 

(11) 

(12) 


(13) 


(14) 


(15) 


at 


h(x,y,x)  -  0  , 

uh  ♦  B(vh  +  wh  )  «  0 


(16) 


53.  DERIVATION  OF  THE  EQUATION  OF  THE  K-dV  TYPE 


-6“ 


■object  to  tha  boundary  conditions 


•“n*  +  kUini5  '  W2  '  W1zni  +  nit  =  °  '  (36) 

P2  -  n2  -  at  z  =  d,  -b  <  y  <  b2  (37) 

n2  -  0  at  z  ■  d,  y  =  -b^b^  (38) 

vji  +  wji  “  -u.h  at  h  «  0  .  (39) 

2T  y  Z  z  lx 

By  differentiating  (34),  (35)  with  respect  to  y  and  z  respectively,  adding,  and  making 
use  of  (19)  and  (20),  we  obtain 

?2P2  “  -*2PiK  •  (40) 

Then  we  differentiate  (22)  and  (25)  with  respect  to  S  and  make  use  of  (34)  and  (35)  to 


obtain 


p2*  -  -“Vise  at  1  “  d'  _bi  <  y  4  b2 ' 


p_  h  +  h.h  ”  0  at  h  «  0 
2y  y  2z  e 


P2  “  -  ♦<t,x,y,z)pl55  +  X2(t,x,S)  . 


(40)  to  (42)  imply 


V2*  -  k2  in  D  , 


“  m2^  at  z  «  d,  -b1  <  y  <  b^  ,  (45) 

♦  h  +  ♦  h  -0  at  h  -  0  .  (46) 

y  y  z  z 

The  Neumann  problem  posed  by  (44)  to  (46)  is  solvable  as  a  consequence  of  (28).  If  we 
assume  ui'vi'wi'P'|'Pig  tend  to  zero  as  S  tends  to  infinity,  it  follows  from  (20),  (34), 
(35)  and  (43)  that 


(“,,»,,*,)  -  w  (kp,'  'Vie-  "♦zpie)  * 


Finally,  we  shall  obtain  a  solution  of  n2  in  terms  of  p,,  and  Aj.  From  (26), 


(37),  (38)  and  (43),  h2  satisfies 


’’zyy  "  <1/t),,2  “  "  Vies  *  (VT,*<Z  “  d>piee  "  (1/T,A2  * 


n2  -  0  at  z  -  d,  y  -  -b1#b 


1'  2  * 


-8- 


(50) 


/  h  (h2  +  h2)  ,/2ds  is  ths  lins  integral  over  the  boundary  of  h  »  0  in  each  cross 
j*  x  y  z 

section  (Figure  1),  and  (h  ,h  ,h  )  is  in  the  outward  normal  direction.  Along  a  ray  we 

x  y  z 

let  a  ■  t  and 

d/d a  -  3/9t  +  (dx/dt)  9/9t  -  9/8t  +  6(x)  9/9t  . 

In  terms  of  a,  (54)  becomes 

Vi®  +  Vi  +  Vipie  +  V««  ’  0  *  (56) 

We  may  also  use  x  as  a  variable  along  a  ray,  then  (56)  assumes  the  form 

-  0  . 


Vlx  *  Vi  +  VlpU  *  Vl«? 


(57) 


{4.  A  RECTANGULAR  CHANNEL  WITH  VARIABLE  WIDTH 

Ha  conaidar  a  ayaMtric  rectangular  ehannal  with  conatant  dapth  and  variable  width. 

Each  croaa  aection  D  of  tha  channel  ia  bounded  by  a  “  0,  a  *  1,  y  ■  ±A(x),  where 

i(x)  ia  a  aaooth  function  aatiafying  0  <  (  A(x)  <  *2  ,  and  are  two 

conatant  a.  The  equationa  for  ♦  ,  (44)  to  (44),  now  be  com 

V2*  -  k2  in  D  , 

♦g  -  *»2wi  at  a  -  1,  -A  <  y  <  A  , 

♦y  -  0  at  y  -  -A,A,  0  <  a  <  1  , 

♦B  ■  0  at  a  »  Of  -A  <  y  <  ft  , 

Let  ♦  •  ♦  +  k2a2/2.  Then  +  aatiafiaa 

V2+  -  0,  In  0  i 

♦a  -  -k2  ♦  w2*,,  at  a  -  1,  -  A  <  y  <  A  , 

♦y  ■  0,  it  y*  -A#A#  0  <  a  <  1  , 

♦B  -  0,  at  8*0,  -A  <  y  <  A  . 

which  nay  be  eolved  by  eaparation  of  varieblea ,  where  by  (27) 

Vj  -  1  -  ooahvy/coahuA  ,  (58) 

and  k,w  are  aaauawd  to  be  given  by  Of). 

By  ay—etry,  we  My  oonaider  the  problen  in  0  <  y  <  A  and  preacribe  the  condition 
♦y  »  0  at  y  *  0,  0  <  a  <  1.  It  ia  aaaily  found  that 

4  -  k2a2/2  *  l  h  coah  nva/A  coa  nwy/A  ,  (59) 

n-0  n 

A  -  (-1)n+1*iw2tanhvA/[nv(u2  ♦  (n«/A)2)ainh  n»/ A),  n  >  1  , 
n 

2 

where  AQ  ia  arbitrary  and  wa  chooae  Afl  -  -  k  /2  ♦  T.  Next  we  proceed  to  find  the 
aolution  for  v  .  Proai  (51),  (52),  (58)  and  (59),  we  have 


1 


2  2  r 

».  -  u  v_  »  coshp  y/coshp 1  +  p  Y  A  cosh  n»/l  cos  n*y/t,  -t  <  y  <  t  , 

2yy  2  n“,  » 

■  0  at  y  -  -1,1  • 


By  syaiMtry  again,  we  may  solve  for  the  half  interval  0  <  y  <  1  and  prescribe 

y„  «  0  at  y  «  0.  Me  find  that 
■*2y  ' 

»  -t  tanh  lit  cosh  py/[2p  cosh  pi)  ♦  y  sinh  py/[2p  cosh  lit] 


+  I  B  ((-1)Bcosh  yy/cosh  ill  -  cos  n*y/t)  . 
.  n 
n  •  1 


(60) 


where 


B  -  p2A  cosh  n»/t  /  III2  +  (nw/1)2]  . 
n  n 

The  expression  for  4  and  given  by  (59)  and  (60)  enable  us  to 

coefficients  in  (54)  for  a  channel  of  constant  depth  but  variable  width, 
we  choose  the  positive  sign  for  6.  Zt  is  obtained  froai  (55)  that 

m0  -  41  [1  -  (1/|il)tanh  pt]  , 

mt  -  21  [1  -  (1/tiDtanh  pt] 1/2  . 

mj  -  ( 1/2)d  a^/dx  , 

mj  -  6kt  [1  -  (1 /lit) tanh  lit] 3/2 


evaluate  the 
To  be  definite. 


♦  J  [4«nw/(ii3t2)]  [1  ♦  (n*/tp)2]”2tanh2pt  tanh  n*/l  , 
n-1 

m4  -  -(2t2k3/3)l1  -  (1/pDtanh  pt] 1/2 

♦  (wt/p2) (tanh2pl  -  1)  ♦  (»/p3)tank  pt 

m 

+  J  (4*»3p2/nw)  [p2  ♦  (n*/l)2]“2tanh  pi 
n«1 

Here  we  keep  k  and  u  in  a^  and  a4  since  other  choices  of  8  are  also  possible.  He 
also  remark  in  passing  that  all  termwise  differentiations  and  integrations  needed  in  this 
derivation  of  m^  to  m^  can  be  justified  without  difficulty,  and  is  independent  of 
the  arbitrary  coefficient  AQ  in  ♦  because  of  (28). 
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$5.  DISCUSSION 

We  observe  f roe  the  Hamilton- Jacobi  equation  (28)  that  the  surface  tension  of  the 
liquid  has  the  apparent  effect  of  reducing  the  width  and  increasing  the  mean  depth  of  the 
channels  Therefore,  the  group  velocity  of  a  capillary  wave  is  always  greatr  than  that  of  a 
surface  wave  without  surface  tension,  and  is  an  increasing  function  of  the  nondiaensional 
surface  tension  coefficient  T.  One  solution  for  (28)  is  given  in  (30).  On  the  other 
hand,  we  may  prescribe  S  -  x  at  t  -  0.  Then  S  is  determined  implicitly  by 

x 

/  10(C) J~dC  -  t. 


In  this  case. 


u  -  G(8),  k  -  G(s)/G(x)  . 


The  disadvantage  of  determining  S  Implicitly  may  be  compensated  by  the  fact  that  the 
initial  condition  for  (34)  at  t  *  0  can  be  directly  expressed  in  terms  of  (  «  '5. 

Certainly  one  can  also  solve  (28)  by  prescribing  data  on  any  simple  smooth  wtve 
6  *  fc0  “  xo  “  in  t,x-plane  where  we  may  identify  S  with  s.  Then 


8  satisfies 


1 

/  (G<n>)dn  -  t  -  tn (s)  . 


Zt  is  shown  that  the  rays  do  not  intersert  each  other  if  G  has  no  characteristic 
direction. 

We  may  derive  the  so-called  Green's  law  (Lamb,  1932)  for  the  change  of  wave  a^litude 
as  follows.  If  we  linearised  the  governing  equations,  and  applied  the  ray  method  (Keller 


(1938),  we  would  obtain  from  (56) 


■Vix  +  Vi  - 0 


along  a  ray.  Therefore, 


P,  *  (PjJqMPI-  /  (i^/m^dc)  , 


where  (p^ )Q  is  the  initial  value  of  p  ,  and  by  (26) 
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(60) 


i 


X 

n,  -  vl(xfy)(p1>0«9(-  /  (■2/*i1)d))  . 

X0 

W*  not*  that  P1,ni  should  b*  lnt*rpr*t*d  r**p*ctly*ly  as  tha  wav*  amplitude  of  tha  first 
ordsr  prasaura  and  surfac*  alavatlon.  Tha  axponantial  in  (60)  will  b*  called  tha 
anplification  factor*  As  scan  fraa  (55),  n^  And  m^  do  not  dapand  upon  tha  solution  0 
of  tha  Neumann  problem  posad  by  (44)  to  (46),  and  can  always  b*  datarminad  for  a  given 
channel.  For  a  rectangular  channel  with  variable  width  b(x)  and  depth  d(x),  it  is  found 
that 

*»,  -  vl(x,y)(p1)(J[b(1  -  (1/ub)tanh  iib/2)]“1/2d"1/4  . 

We  see  that  the  amplification  factor  is  always  greater  for  a  surface  wave  with  surface 
tension. 

If  the  channel  has  a  uniform  cross  section,  then  h^  -  (G_1 )x  ■  0  and  from  (55) 

“  0,  and  all  other  coefficients  are  constant.  Furthermore,  the  rays  are  straight 
lines  and  given  by  8  -  t  ■  constant.  We  obtain  from  (50)  that 

Vic  ♦  Vl>«  ♦  View  “  0  * 

from  which  expressions  for  progressive  waves  of  permanent  type  can  be  found.  Finally  if 
T  tends  to  aero,  that  is,  V  tends  to  infinity,  (55)  reduce  to  those  obtained  in  Zhong 
and  Shen  (1982)  for  waves  without  surface  tension. 
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